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Supersonic cavity flows are characterized by compression and expansion waves, shear layer, and oscillations in- 


side the cavity. For decades, investigations into cavity flows have been conducted, mostly with flows at zero 


pressure gradient entering the cavity in straight walls. Since cavity flows on curved walls exert centrifugal force, 


the features of these flows are likely to differ from those of straight wall flows. The aim of the present work is to 


study the flow physics of a cavity that is cut out on a curved wall. Steady and unsteady numerical simulations 


were carried out for supersonic flow through curved channels over the cavity with L/H = 1. A straight channel 


flow was also analyzed which serves as the base model. The velocity gradient along the width of the channel was 


observed to increase with increasing the channel curvature for both concave and convex channels. The pressure 


on the cavity floor increases with the increase in channel curvature for concave channels and decreases for con- 


vex channels. Moreover, unsteady flow characteristics are more dependent on channel curvature under supersonic 


free stream conditions. 
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Introduction 


The subject of transonic and supersonic flows over a 
wall mounted cavity has long been of research interest 
for many engineers and scientists. Although a wall moun- 
ted cavity is geometrically very simple, it involves many 
diverse flow characteristics such as shear layer, shock 
wave, flow separation, recirculation regions, etc. Much 
has been learned from previous computational as well as 
experimental research [1-3]. For a given cavity geometry, 
both the Mach number as well as the Reynolds number 
upstream of the cavity have important influence on the 
flow field inside the cavity. Depending on the shape of 
the cavity as well as its length to depth ratio (L/H), flow 
over a cavity can be specified as an open cavity flow, a 
transitional flow or a closed cavity flow. For such cases, 
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the cavity drag and flow oscillation characteristics are 
well documented [4]. 

With the advent of advanced computing setups, efforts 
have been made to study the oscillatory flow mechanisms 
employing various numerical simulation techniques [5-6]. 
Flow oscillations within the cavity occur predominantly 
in the spanwise, longitudinal or transverse directions de- 
pending on the length to depth ratio of the cavity. Ac- 
cording to Plumblee et al. [7] the turbulent shear layer 
which spans the free boundary of the cavity provides the 
cavity with a broad band noise source from which the 
cavity detects certain frequencies to amplify. This me- 
chanism is verified for cavities exhibiting transverse mod- 
es of oscillation which exist when the cavity length to 
depth ratio is equal to or less than 1. As the cavity be- 
comes shallow (L/H > 1), the oscillation mechanism 
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Nomenclature 
A surface area (m°) U 
Cp cavity drag (N) u 
Cp pressure coefficient xX, y 
D drag (N) y 
f frequency (Hz) 
H cavity height (mm) T 


ratio of convective velocity of vor- 
tices to the freestream velocity 
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free stream velocity (m/s) 
axial velocity (m/s) 
Cartesian coordinates 


nondimensional wall distance 


Greek letters 


shear stress 


a factor to account for time lag between the passage 
of a vortex and emission of acoustic pulses at the 
trailing edge 


stagnation condition 
center of the cavity floor 
cavity floor 

inlet condition 

leading edge 

static state 


root mean square 


L cavity length (mm) Subscripts 
M Mach number 0 

m mode number (1,2,3,...) c 

Ẹ pressure (Pa) floor 

R radius of curvature (mm) in 

SPL sound pressure level (dB) leading 
Stn Strouhal number st 

T temperature (K) rms 

t time (s) trailing 


shifts from an acoustic response of the cavity to an acou- 
stic feedback mechanism [4]. Extensive research and 
reports have been presented on flow over cavities exhi- 
biting acoustic feedback mechanism as these types of 
cavities are practically more applicable [8]. Investiga- 
tions have pointed out that for cavities of L/H > 1, the 
driving mechanism of oscillation is the property of flow 
over the cavity rather than being the property of the ap- 
proaching boundary layer. However, for deep cavities 
(L/H < 1), the driving mechanism is the shear layer ef- 
fects which develop from the incoming boundary layer. 
Heller and Bliss [9] observed that the transition of the 
incoming boundary layer from a laminar to a turbulent 
one more significantly affects the cavity response than 
that for a turbulent upstream boundary layer. The pres- 
sure fluctuations in the cavity are insensitive to the change 
in the Reynolds number. Gharib and Roshko [10] ob- 
served a significant change in the behavior of cavity os- 
cillations when the ratio of cavity length relative to the 
upstream boundary layer momentum thickness was in- 
creased during their incompressible experiment for an 
axisymmetric cavity. Investigations conducted by several 
investigators [11-13] underlined the effects of the in- 
coming boundary layer on the shear layer growth over 
the cavity. The incoming boundary layer separates at the 
leading edge of the cavity. It then transforms into a free 
shear layer and interacts with the trailing edge depending 
on various geometric and fluid dynamic aspects. The 
shear layer behavior governs the different flow features 


trailing edge 


exhibited in the vicinity of the cavity. 

Most of the works completed up to now have focused 
on transonic and supersonic flow over the cavities moun- 
ted on a straight wall where the ratio of the radius of the 
wall curvature to the cavity depth is infinite; the incom- 
ing flow has typical laminar or turbulent boundary veloc- 
ity profiles. However, many engineering and industrial 
applications have been made with cavities over curved 
wall surfaces such as vortex combustors with cavities, 
external cutouts on aircraft bodies, pipeline-cavity sys- 
tems, etc. In high-speed flows over a curved surface, the 
incoming flow upstream of the cavity is subjected to a 
centrifugal force with its strength depending on the ra- 
dius of wall curvature and Mach number. The centrifugal 
force induces a pressure gradient opposite to its direction, 
leading to a variation in the velocity profile of the up- 
stream flow. Under such conditions, the characteristics of 
the cavity flow will considerably differ from those pre- 
sented in many of the previous works involving a straight 
wall. To the best of the present authors’ knowledge, no 
investigation has been made on high-speed flows over a 
curved wall mounted cavity. In the present study, a com- 
putational fluid dynamics method was applied to investi- 
gate high-speed flows over rectangular cavities with an 
L/H value of 1. The cavities are installed in curved chan- 
nels by varying the radius of curvature, with the L/R ra- 
tios from 0 to 0.11. Both concave and convex channels 
are used for the present study. This facilitates the study of 
flow characteristics in the channel under both the positive 
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and negative effects of the centrifugal force. The setup is 
such that the cavity mounted on a concave channel is 
subjected to a positive centrifugal force, whereas that on 
a convex channel undergoes a negative centrifugal force. 
The Mach number in the channel inlet varies from 0.4 to 
1.8, with the Reynolds number based on the channel 
height varying from 0.56 x 10° to 3.21 x 10°. 

The present computational fluid dynamics (CFD) re- 
sults are validated with the existing experimental data. It 
is found that the cavity pressure loadings and cavity drag 
depend heavily on the channel curvature. It is observed 
that the unsteady pressure loadings depend on the type of 
curvature when compared with the base model. 


Numerical Analysis 


Test model and numerical methodology 

The test domain consists of a 58.42 mm wide channel 
with a rectangular cavity mounted on the wall. A straight 
walled channel as shown in Fig.1 (a) is used as the base 
model and six other channels with curved walls as shown 
in Figs. 1(b) and 1(c) are used for the study. The width of 
the channel is fixed through an iterative procedure as it is 
observed that the channel width plays an important role 
in establishing a uniform flowfield in the vicinity of the 
inlet region for the curved channels. For a width larger 
than 58.42 mm, the pressure conditions imposed at the 
channel exit tend to slow down the flow at lower Mach 
numbers. For a smaller width, the flow tends to become 
stagnated inside the channel at higher Mach numbers. 
The channel width is thus fixed such that the freestream 
flow properties remain similar for the various models at 
different inlet conditions. The cavity is located approx- 
imately 100 mm downstream of the inlet section of the 
channel for all the models. The dimensions of the cavity 
are L = 25.4 mm and of H = 25.4 mm and it therefore has 
an L/H ratio of 1. The flow channel curvatures are 
changed by altering the radius of the curve as shown in 
Fig. 1, which thereby provides three L/R ratios of 0.04, 
0.08 and 0.11. The straight channel is considered to have 
an infinite radius and hence an L/R ratio of 0.0 and no 
centrifugal effects. The geometrical aspects of the differ- 
ent models are provided in Table 1. A uniform flow field 
at Min = 0.4, 0.8, 1.4 and 1.8 is established at the inlet of 
the flow domain. 

The analysis is performed using the commercially 
available computational fluid dynamics package, An- 
sys15. The two-dimensional (2D) flow domain is gener- 
ated using the modeling software, ICEM-CFD. Finer 
meshes are used adjacent to the walls in order to solve 
the boundary layer interactions. 1300 x 383 mesh ele- 
ments are used for the base model, which gives a wall 
y ~1 for the converged steady solutions at a location 5 
mm upstream of the cavity on the bottom wall of the 


flow channel. Both steady and unsteady Reynolds Aver- 
aged Navier-Stokes equations are used to simulate the 
flow field in the channel. As the flow over the cavity is 
presumed to be predominantly oscillatory, the unsteady 
scheme is applied on the converged steady solutions with 
a time-step of At = 10° s. Based on the theoretical pre- 
dictions by Heller et al. [9], this time-step is deemed to 
be sufficient for resolving the flow periodicity in the cur- 
rent model for the flow at different M;„. Pressure inlet and 
outlet boundary conditions are used at the entry and exit 
of the flow domain respectively. Other domain bounda- 
ries are designated as adiabatic walls with no-slip boun- 
dary conditions. The inlet total temperature is kept at To= 
300 K for all the test cases. Multiple viscous turbulence 
schemes such as k-s and k-w and their variations are used 
to model turbulence, and the flow analysis results are 
compared with the available experimental data [14] for 
the purpose of validation using the base model. Fig. 2 
shows the variation of the coefficient of pressure on the 
cavity floor for the different turbulence models used in 
the flow analysis of the base model at M;,= 1.5. It is ob- 
served that the k-æ SST model offers the closest proxim- 


(a) Schematic of base model 


LEJ 


(b) Concave channel 


(c) Convex channel 
Fig.1 Schematic diagrams of test setup 


Table 1 Details of the test models 


L R L/R Model type 
0 Base model 
0.04 Concave/Convex 
25.4 œ - 222.6 
08 Concave/Convex 
0.11 Concave/Convex 
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ity to the experimental results. Numerical simulations are 
carried out with the k- turbulence model for the remain- 
ning test cases. To check the dependence of flow para- 
meters on the flow domain mesh, a grid independence 
study is carried out with three different grids using the 
base model at M;,= 1.8. The grid is clustered adjacent to 
the walls to fully resolve the boundary layer effects. The 
time averaged static pressure at the cavity floor for the 
different grids used in the study is shown in Fig. 3. The 
static pressure is observed to vary as the grid inside the 
cavity is changed from coarse (120 x 94) to fine (200 x 
156). The grid in the main flow channel is also altered in 
proportion with the cavity grid. However, as the compu- 
tational time for the fine grid is much longer than that of 
the coarse grid, the intermediate grid of 160 x 125 ele- 
ments inside the cavity is used. For this grid, the static 
pressure distribution at the cavity floor is almost equal to 
that obtained during computation with the fine grid. The 
periodicity of flow oscillations are determined from the 
unsteady analysis of the converged steady flow solutions 
for which the mass imbalance between the inlet and out- 
let falls below the order of 10°. 


Leading edge Cavity floor Trailing edge 
N GŘ 


© Exp. Data [Ref.14] 
—k-epsilon standard 
sasa k-epsilon RNG 
— —k-epsilon realizable 
— --k-omega standard 
— -- k-omega SST 


Fig. 2 Comparison between experimental results and CFD 
prediction (L/H = 1 and M,,, = 1.5) 
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Fig. 3 Static pressure variations at cavity floor with different 
grids 
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Results and Discussion 


In order to determine the incoming flow field and the 
boundary layer structure, the average of the velocity gra- 
dient values along the width of the channel at a location 
approximately 28 mm upstream of the cavity leading 
edge is plotted for different models as shown in Fig. 4. 
For the convex models, the flow velocity increases along 
the floor of the channel and decreases along the top wall 
of the channel as compared to the base model. The con- 
cave channels showed a reverse trend compared to that of 
the convex channels. The maximum flow velocity is ob- 
served to occur near the top wall of the concave channel 
with the minimum velocity occurring near the bottom 
wall. If the channel is considered to be a part of a much 
larger circular flow field, the incoming velocity distribu- 
tion gives an approximation of a vortical flow that occurs 
over the cavity. 

To differentiate between the concave and the convex 


models, the (@u/@y) values were considered to be pos- 


itive for the concave channels and negative for the con- 
vex channels. The velocity distribution is obtained from 
the time averaged steady results of the numerical simula- 
tion. For a uniform flow of M= 0.4, 0.8, 1.4 and 1.8 at 
the inlet of the channel, the velocity distribution along 
the channel width shows dependency on the curvature of 
the channel as shown in Fig. 4. As the curvature increas- 
es, which is indicated by an increase in the L/R ratio, the 


(Gu / dy) value increases for the same M,,. As a result, 


the shearing effects on the flow field will increase at a 
higher channel curvature. 
The average velocity gradient is also observed to in- 


crease as the centrifugal effects on the flow field increase. 


An increase in the flow of Min increases the average ve- 
locity gradient. At higher L/R ratios, the effect of change 
in M; is more pronounced than at the lower L/R ratio. It 


5000 


Closed symbols : Concave channel 
Open symbols : Convex channel 

—+—M,,=1.8 
-4 M,,=0.8 


4000 


3000 = 
-& M,=0.4 


2000 
1000 


cu 


0 0.04 0.08 0.12 


L/R 


Fig.4 Average velocity gradient upstream of the cavity (x/L = 
1.1) 
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is also observed that as the curvature of the flow channel 
increases from L/R = 0.04 to L/R = 0.11, the non- 
uniformity of the velocity field also increases. 

Since the incoming flow field over the cavity changes 
due to the channel curvature, it is important to estimate 
the pressure loadings on the cavity. Static pressure load- 
ing on the cavity and its dependency on the flow path 
curvature is investigated by plotting the time averaged 
static pressure at the center of the cavity floor for all the 
test models at different M. To compare the static pres- 
sure for different models, the static pressure has been 
non-dimensionalized using the static pressure at the 
channel inlet. The plots are shown in Fig. 5. 
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Fig.5 Time averaged static pressure at the centre of the cavity 
floor 


The P, loading in non-dimesionalized form at the floor 
of the cavity for the base model shows little variation for 
different Mm, as shown in Fig. 5. Quantitatively, the 
loading is 101 kPa for M;, = 0.4 and 99 kPa for Min =1.8. 
As the curvature of the concave channels increase, the P. 
increases. The increase is higher at supersonic flow than 
for the subsonic flow cases. The P, for the L/R = 0.11 
concave channel at M,,=1.8 is almost twice that for the 
L/R = 0.0 channel. For the convex channels, the P, on the 
cavity floor is observed to decrease with an increase in 
Min and channel curvature. The pressure loadings on the 
floor of the cavity mounted on the convex channel wall 
are found to be lower than those of the base model. A 
minimum pressure loading is observed for the Z/R = 0.11 
convex channel among all the models. Earlier studies 
indicate that, as the fluid mass exchange between the 
cavity and the freestream increases, higher levels of 
pressure loadings are observed within the cavity [15, 16]. 
Consequently, the fluid mass exchange rate between the 
cavity and the freestream is investigated for all models. 
The time averaged mass exchange between the cavity 
and the mainstream flow is calculated by integrating the 
mass flow rate over an imaginary line joining the leading 


edge and the trailing edge of the cavity which is parallel 
to the channel walls. It is observed that at M;, =1.8 for the 
L/R = 0.11 concave model, the mass flow rate of 55.47 x 
10° kg/s is the maximum among all the models analyzed. 
Conversely, the mass flow rate of 10.04 x 10° kg/s for 
the L/R = 0.11 convex model is observed to be the mini- 
mum among all the models. The time averaged static 
pressure level at the cavity floor as shown in Fig. 5 seems 
to depend on the fluid mass interaction between the cav- 
ity and the freestream. It is also observed that the varia- 
tions of P, for the different models are similar to those 
shown in Fig. 4. The increase in the value of the average 
velocity gradient indicates greater shearing effects in the 
flow. This amplification of shearing effects may alter the 
boundary layer and hence the transport of fluid mass and 
momentum into the cavity, contributing to higher levels 
of pressure at the cavity floor. 

As the basic fluid dynamic configuration being inves- 
tigated is also applicable to the external cutouts of air- 
crafts, investigating the drag generated due to the pres- 
ence of such configurations is an integral part of such 
studies. Considering internal flow and applicability to 
vortex combustors, the drag generated in the flow field 
provides an estimation in relation to the loss of flow 
momentum. Fig. 6 shows a comparison between the Cp 
generated in the flow channels of varying curvatures at 
different flows of M;, and the Cp generated for the base 
model at the respective Min. The Cp is calculated based 
on equation (1) [17]. 

D= | PdA - È PdA + | a4 (1) 
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Fig. 6 Variation of cavity drag for the different models (a = 
concave +, convex -) 


Here, the channel curvature is differentiated between 
concave and convex by representing the L/R ratio as pos- 
itive for the concave channels and as negative for the 
convex channels. For both the concave and the convex 


150 


models, the Cp generated is observed to increase with an 
increase in channel curvature at all M;,. From Fig. 6, it is 
observed that as M;, increases, the ratio of Cp/Cpy/r=0.0) 
also increases for the same model. This increase is more 
prominent at supersonic Min. The overall drag increase 
with channel curvature could be an effect of the higher 
shearing force depending on the average velocity gra- 
dient for the respective channel. 

The unsteady pressure variations at the centre of the 
cavity floor are shown in a transformed power spectral 
density (PSD) vs frequency domain in Figs. 7(a) and 7(b) 
for the different models at M;,=0.8. As the cavity L/H 
ratio is 1, the unsteady pressure is measured at the centre 
of the cavity floor, since transverse oscillations are pre- 
dominant in these types of configurations [7]. For the 
base model, the dominant oscillation mode is observed to 
be at 1.038 kHz at M,,,= 0.8. For the concave models, the 
frequency of the dominant mode is observed to remain 
almost the same as that of the base model. The frequency 
of the dominant mode is observed as 1.099 kHz only for 
the Z/R=0.11 concave model. However, the sound pres- 
sure level (SPL) of 20 log (Prms/20 Pa) (reference pres- 
sure of 20 uPa) of the dominant mode decreases as the 
curvature of the channel increases for the concave mod- 
els. For the convex channels, at M;,= 0.8, the frequency 
of the dominant mode is also observed to remain almost 
the same for all models. The frequency of the dominant 
mode for the L/R = 0.04 convex model is observed to be 


PSD 


10 
8 3 
LIR 0.04 4 6 x 10° 
2 
(Base model)0 0 ~ Frequency(Hz) 


(a) Concave channel 
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Fig.7 PSD of the fluctuating pressure at the centre of the cav- 


ity floor (Mi, = 0.8) 
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the same as that of the base model at 1.038 kHz. For the 
L/R = 0.08 and L/R = 0.11 convex models, the frequency 
of the dominant mode is observed to be 1.099 kHz. 
Hence, no remarkable change is observed in the domi- 
nant mode frequency for the different channels with va- 
rying curvatures at M;,= 0.8. The magnitudes of the do- 
minant modes in terms of SPL for the convex channels 
are observed to be almost the same as that of the base 
model. The oscillation frequencies at M;,= 1.8 are ob- 
served to be very high for all the models. The pressure 
spectra for the different models are shown in Fig. 8. 

For the L/H = | cavities, several previous works indi- 
cate that the resonant frequencies are much higher than 
the cavities with L/H > 1 with supersonic flow. The 
pressure spectra for the concave models are marked 
mostly by broadband fluctuations for lower L/R ratio 
channels. As the curvature of the concave channels in- 
creases, the discrete modes are observed. The convex 
models oscillate at discrete peaks rather than broadband 
oscillations as the channel curvature increases. For the 
baseline cavity of Z/R = 0.0, the first oscillation mode 
occurs at 47.2 kHz and the dominant mode occurs at 92.1 
kHz. These modes correspond to mode 5 and mode 10 of 
the theoretically calculated mode frequencies, respec- 
tively [9]. For the L/R = 0.04 concave model, apart from 
the broadband fluctuations, three other discrete modes 
are observed and the dominant mode corresponds to a 
frequency of 85.45 kHz. This is within 100 Hz of the 


(a) Concave channel 


(b) Convex channel 


Fig. 8 PSD of the fluctuating pressure at the centre of the 


cavity floor (Mi, = 1.8) 
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theoretically calculated frequency of mode number 9. As 
the curvature of the concave channels increases, the do- 
minant mode shifts to a lower frequency with the sup- 
pression of pressure fluctuations. The dominant ampli- 
tude of oscillation for the L/R = 0.0 model is observed to 
be higher than the amplitude of the curved concave 
channel models. 

For the convex models, high amplitude discrete modes 
are observed along with broadband fluctuations. The 
pressure spectra of the convex models are shown in Fig. 
8(b). As the channel curvature increases, the oscillations 
move from broadband type to discrete modes. As the 
curvature of the flow channel increases, further high am- 
plitude discrete peaks are observed in the pressure spec- 
tra. The dominant modes for the Z/R = 0.08 and L/R = 
0.11 convex models are observed to be at a frequency of 
63.78 kHz and 227.2 kHz, respectively. The sound pres- 
sure levels for the convex channels are observed to be 
higher than those for the concave channels, and this in- 
dicates a suppression of oscillation for the concave 
channels using supersonic flow. The SPL of oscillations 
for M= 1.8 flow is much lower than that for the M, = 
0.8 flow for the curved channels. 

In terms of quantitative values of the amplitude of os- 
cillations, it is observed that for subsonic Min = 0.8, the 
values are similar to those for the base model. Previous- 
ly, (6u/ éy) at subsonic M; was observed to vary slightly 
between the base model and the curved channel models, 
Ou / Oy 


as shown in Fig. 4. At supersonic Min, is high- 


er and the amplitudes of pressure oscillations are subdued. 


Higher average velocity gradients and shearing effects 
are thus observed to attenuate the oscillation amplitude in 
the flow field. 

The mode frequencies at the two M; of 0.8 and 1.8 are 
compared with the theoretically predicted mode frequen- 
cies using Rossiter’s modified formula [9]. The modified 
Rossiter’s formula is given as 

JL m-E 
Stn = 7 = 1, 7 i (2) 
[l+(y-)M2/2] * 

The values of the empirical constants € and 1/K are 
evaluated through the best fit of the measured data as ¢ = 
0.25 and 1/K = 1.75. The theoretically and numerically 
obtained Strouhal numbers (Stn) are plotted for different 
Min and are shown in Fig.9. For the subsonic M= 0.8, 
the oscillation frequencies are observed to be concen- 
trated in the first four fundamental modes. At supersonic 
Mi, = 1.8, Stn are mostly spread over high frequency 
modes. It is observed from Fig. 9 that for the curved 
channels, the oscillation frequencies do not agree with 
the predicted values at higher M. Quantitatively, the 
oscillation frequencies exhibited by the base model agree 
more closely with the predicted values. For the curved 


channels, at supersonic M; the variations from the pre- 
dicted frequencies are observed to occur for all the modes. 
At Mn= 0.8, the oscillation frequencies do not vary con- 
siderably with the channel curvature. However, as the 
channel curvature increases, the departure from the pre- 
dicted values of oscillations is observed to increase. For 
the lower modes, the frequencies are observed to be 
within 50 Hz of the predicted values. The shift in oscilla- 
tion frequencies from the predicted values is higher at 
higher modes. 


—@Concave 
-> Convex 
A Base model 
---- Heller & Bliss[9] 


—@ Concave 
=- Convex 
A Base model 


(b) M, = 1.8 


Fig.9 Variation of Stn for different models 
Conclusions 


Steady and unsteady numerical simulations were car- 
ried out for the flow through curved channels over an 
L/H = 1 cavity for both subsonic and supersonic flows. 
Flow through a straight channel was also analyzed and 
served as the base model. Time averaged pressure load- 
ings on the cavity are observed to increase with an in- 
crease in the channel curvature for concave channels and 
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a decrease in the channel curvature for convex channels. 
As a result of the fluid mass interaction between the cav- 
ity and the freestream, the pressure loadings are higher 
than those for the base model when the cavity is sub- 
jected to positive centrifugal force, and are lower when it 
is subjected to negative centrifugal force. The effect of 
channel curvature on the steady pressure characteristics is 
greater at a higher Mach number than at lower Mach 
number. The quantitatively larger fluid mass exchange 
between the cavity and the mainstream flow at a higher 
Mach number seems to be the driving mechanism behind 
this observation. The drag in the flow channel increases 
at higher L/R ratio for both concave and convex models. 
The drag for the curved channels is higher than that for 
the straight channel. The flow field is observed to be 
dominated mostly by broadband oscillations for the con- 
cave models and discrete oscillations for the convex 
models at a supersonic Mach number. Unsteady flow 
characteristics are observed to be more dependent on 
channel curvature at supersonic Mach numbers. A shift in 
the frequency of the dominant mode is observed with a 
change in channel curvature. The suppression of oscilla- 
tory flow at the supersonic Mach number is observed for 
the concave channels. Frequencies of oscillation modes 
of the curved channels show a departure from the pre- 
dicted values of the straight channel. 
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